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Multiferroic devices hold profound promise for ultra-low energy computing in beyond Moore’s law era. The magnetization of a 
magnetostrictive shape-anisotropic single-domain nanomagnet strain-coupled with a piezoelectric layer in a multiferroic composite 
structure can be switched between its two stable states (separated by an energy barrier) with a tiny amount of voltage via converse 
magnetoelectric effect. With appropriate choice of materials, the magnetization can be switched with a few tens of millivolts of 
voltages in sub-nanosecond switching delay while spending a miniscule amount of energy of ~1 attojoule at room-temperature. 
Here, we analyze the area-delay-energy trade-offs of these multiferroic devices by solving stochastic Landau-Lifshitz-Gilbert equation 
in the presence of room-temperature thermal fluctuations. We particularly put attention on scaling down the lateral area of the 
magnetostrictive nanomagnet that can increase the device density on a chip. We show that the vertical thickness of the nanomagnet 
can be increased while scaling down the lateral area and keeping the assumption of single-domain limit valid. This has important 
consequence since it helps to some extent preventing the deterioration of the induced stress-anisotropy energy in the magnetostrictive 
nanomagnet, which is proportional to the nanomagnet’s volume. The results show that if we scale down the lateral area, the switching 
delay increases while energy dissipation decreases. Avenues available to decrease the switching delay while still reducing the energy 
dissipation are discussed. 


Index Terms —Nanoelectronics, spintronics, multiferroics, energy-efficient computing, straintronics, area-delay-energy trade-offs. 


I. Introduction 

S TRAIN-MEDIATED multiferroic devices, i.e., a magne¬ 
tostrictive layer strain-coupled to a piezoelectric layer 
works according to the principles of converse magnetoelectric 
effect 0, 0, 0, 0, 0. When a voltage is applied across 
such heterostructure, the piezoelectric layer gets strained and 
the strain is transferred elastically to the magnetostrictive 
layer and generates a stress-anisotropy (or magnetoelastic 
anisotropy) in it. If we consider a shape-anisotropic single¬ 
domain nanomagnet having two stable magnetization states 
as the magnetostrictive layer, the generated stress-anisotropy 
can overcome the shape-anisotropy energy barrier and rotate 
the magnetization. (See Fig. ffl) With appropriate choice of 
materials, the magnetization can be switched between its two 
stable states separated by an energy barrier in sub-nanosecond 
switching delay while expending miniscule amount of energy 
of -1 attojoule at room-temperature 0, 0- Such study 
has opened up a new field called straintronics Q, 0, 0, 
which can possibly replace the conventional charge-based 
electronics as our future information processing paradigm. 
Experimental efforts demonstrating such electric field induced 
strain-mediated magnetization rotation are emerging too 0, 
na, hd, na, in. n should be emphasized that remedies 
for substrate clamping effect particularly for low-thickness 
piezoelectric layers (< 100 nm) need to be considered m, 
nn. The use of thin films rather than thick substrates m 
for piezoelectric layers allows us to work with lower voltages 
and therefore it decreases the energy dissipation. 

The stress-anisotropy induced in the magnetostrictive nano- 
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Fig. 1. A voltage-controlled strain-mediated multiferroic composite device. 
The magnetization of the single-domain magnetostrictive nanomagnet can 
be switched between its two stable states by the generated stress on the 
nanomagnet via applying a voltage across the composite structure. The voltage 
generates strain in the piezoelectric layer, which is transferred elastically to 
the magnetostrictive layer. 


magnet is proportional to the magnetostrictive coefficient of 
the material used, the generated stress (which is in turn 
proportional to the electric field), and the volume of the nano¬ 
magnet. For a given material with a certain magnetostrictive 
coefficient and a fixed volume, a higher stress facilitates having 
a faster switching speed (for a sufficiently fast ramp rate of 
stress) but causes a higher energy dissipation (6). The energy 
dissipation in these multiferroic devices comprises of two 
components: one is due to magnetization damping through 
which magnetization relaxes to a minimum energy position 
dissipating energy and the energy dissipation in the external 
circuitry due to applied voltage (oc CV 2 , where C is the 
capacitance of the piezoelectric layer and V is the applied 
voltage) (3, (6). Since the voltage required to switch the 
magnetization is miniscule, the “CV 2 ” energy dissipation is 
also miniscule. 

In this paper, we study the area-delay-energy trade-offs 
of such multiferroic magnetoelectric devices. We particularly 
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analyze the effect of reducing lateral area to scale down the 
area consumption on a chip. Scaling down lateral area for a 
constant thickness reduces both the shape-anisotropic energy 
barrier height separating the two stable magnetization states 
and the generated stress-anisotropy in the nanomagnet since 
they are both proportional to the volume of the nanomagnet. 
Apparently, the concern is how to deal with such scenario if 
we wish to decrease the lateral area of such devices. It turns 
out that while we reduce the lateral area, we can increase the 
vertical thickness of the nanomagnet to some extent keeping 
the single-domain limit assumption for the nanomagnet intact. 
Increasing thickness while reducing lateral area of the nano¬ 
magnet facilitates preventing a drastic reduction in volume 
and hence the induced stress-anisotropy in the nanomagnet 
does not deteriorate much. Note that we need to keep the 
shape-anisotropic energy barrier constant since that determines 
the error-probability of spontaneous reversal of magnetization 
between two stable states. The decrease in shape-anisotropic 
energy due to decrease in volume needs to be compensated by 
properly adjusting the lateral dimensions. We perform simula¬ 
tions using stochastic Landau-Lifshitz-Gilbert (LLG) equation 
of magnetization dynamics lH71 . (YU, lH9l to extract the 
switching delay and energy dissipation during magnetization 
reversal. 

If we increase the thickness of the nanomagnet, the thick¬ 
ness of the piezoelectric layer needs to be also proportionately 
increased since the effectiveness of strain-transfer depends on 
the ratio of the two layers 0. Modeling the piezoelectric 
layer as a parallel-plate capacitor and noting that the area, 
A of the capacitor decreases while the thickness d increases, 
the capacitance of the piezoelectric layer (oc A/d) decreases 
significantly. The lateral area decreases more compared to the 
increase in thickness of the piezoelectric layer to enforce the 
single-domain limit, i.e., the factor Ad eventually decreases. 
Thus even if we increase the voltage V for a thicker device 
to keep a constant electric field (and thus a constant stress) 
across the device, the “CV 2 ” energy dissipation eventually 
decreases. 

The rest of the paper is organized as follows. Section HU de¬ 
scribes the model utilized to perform the simulations. Stochas¬ 
tic Landau-Lifshitz-Gilbert (LLG) equation in the presence of 
room-temperature thermal fluctuations is solved to trace the 
magnetization dynamics and calculate the associated perfor¬ 
mance metrics, i.e., switching delay and energy dissipation. 
Simulation results showing how the performance metrics area, 
switching delay, and energy dissipation vary with the decrease 
in lateral dimensions (and the increase in thickness) of the 
nanomagnet are presented in Section [Till Finally, Section [TV] 
discusses on the results presented and concludes the paper. 

II. Model 

We model the piezoelectric and magnetostrictive layers as 
elliptical cylinders as shown in Fig. [T] The magnetostrictive 
nanomagnet lies on the y-z plane; the major axis is aligned 
along the ^-direction and the minor axis along the ^/-direction. 
The dimensions of the major axis, the minor axis, and the 
thickness are a, b , and l, respectively. So the nanomagnet’s 


aspect ratio is a/b , cross-sectional area is A = (7r/4)a6, and 
volume is ft = (tt/A) abl. The z- axis is the easy axis, the 
y -axis is the in-plane hard axis, and the x-axis is the out-of¬ 
plane hard axis. Since l < b , the out-of-plane hard axis is 
much harder than the in-plane hard axis. As l becomes higher 
and the lateral dimensions a and b become smaller, it is easier 
for magnetization to deflect out-of-plane. In standard spherical 
coordinate system, 0 is the polar angle and (j) is the azimuthal 
angle of the magnetization vector. Note that when <j) = ±90°, 
the magnetization vector lies on the plane of the nanomagnet. 
Any deviation from (j> = ±90° corresponds to out-of-plane 
excursion. 

We can write the total energy of the magneto strictive 
polycrystalline single-domain nanomagnet when it is subjected 
to uniaxial stress along the easy axis (major axis of the ellipse) 
as the sum of the shape-anisotropy energy and the stress- 
anisotropy energy urn 0, a as 


E = B{4>) sin 2 0 ( 1 ) 

where 

B{4>) = Bshaped 0) T~ Bstress , (2a) 

Bshape (0) = (mo/2 )MgQ[(N d - yy - N d - zz ) (2b) 

4“ {Nd—XX Nd—yy) COS (/>], (2c) 

B stress = (3/2)A s <t£ 2, (2d) 


M s is the saturation magnetization, N d - mrn is the com¬ 
ponent of demagnetization factor along ra-direction, which 
depends on the nanomagnet’s dimensions Go), GD, (3/2)A s 
is the magnetostrictive coefficient of the single-domain mag- 
netostrictive nanomagnet lf20l . and a is the stress on the 
nanomagnet. Note that the potential landscape of the mag¬ 
netostrictive nanomagnet is symmetric in space, however, the 
out-of-plane excursion of magnetization provides an equivalent 
asymmetry to facilitate switching of magnetization towards its 
correct direction when it reaches the x-y plane (9 = 90°) ll22l . 
Thermal fluctuations create an wide distribution of magne¬ 
tization while switching starts and also the time required 
to reach 0 = 90° is a distribution due to thermal fluctua¬ 
tions l22t . Therefore, it necessitates a sensing methodology 
to detect when magnetization reaches around 0 = 90° but 
the sensing does not need to be very accurate as intrinsic 
dynamics provides some tolerance f22l . Note that the product 
of magnetostrictive coefficient and stress needs to be negative 
in sign for stress-anisotropy to overcome the shape-anisotropy. 

From ©, the in-plane barrier height (i.e., when 0 = 90°, 
4> = ±90°, and a = 0) between two magnetization stable 
states ( 0 = 0° and 180°) can be written as 

Ebarrier = (/W 2 )M^Q l (N d - yy - N d - zz ). (3) 

The probability of spontaneous reversal of magnetization 
due to thermal fluctuations is exp[— Ebarrier / kT] according 
to the Boltzmann distribution, where k is the Boltzmann 
constant, and T is temperature. The minimum stress required 
to overcome the barrier height Ebarrier can be determined as 

_ (mo/2 )M 2 (N d . yy - N d . zz ) 
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The minimum stress cr m ^ n does not directly depend on 
volume Cl since both shape-anisotropy and stress-anisotropy 
energies are proportional to volume, however, cr m ^ n depends 
on the particular dimensions of the nanomagnet (i.e., a, b , and 
l) through the dependence of demagnetization factors N d - yy 
and Nd_ zz . It should be emphasized that the barrier height 
E barrier depends on volume Cl and to satisfy single-domain 
limit El) . there are certain dimensions of the nanomagnet that 
we can choose. 

The magnetization M of the nanomagnet has a constant 
magnitude but a variable direction, so that we can represent 
it by a vector of unit norm n m = M/|M| = e r where e r is 
the unit vector in the radial direction in spherical coordinate 
system represented by (r,#,</>). The torque, Te acting on the 
magnetization due to shape and stress-anisotropy is derived 
from the derivative of potential energy as in © 12, 0. There 
is an additional torque, Tth due to room-temperature (300 
K) thermal fluctuations 0- 

The magnetization dynamics under the two aforesaid 
torques is described by the stochastic Landau-Lifshitz-Gilbert 
(LLG) equation El, ED, El as 

dn m f dn m \ I7I r 

__ _ a ^ nm x __j = __ [Te + Tth ] (5) 

where a is the phenomenological damping parameter, 7 is the 
gyromagnetic ratio for electrons, and M = poM s Cl. 

After solving the LLG equation, we get the following 
coupled equations for the dynamics of 6 and f 0 

(l + a 2 ) — = jj:[B shape, <t>((f>)sin0 - 2aB((j))sinO cosO 

+ (ccPq + P(f>)\, (6) 

(l + a 2 ) — = ^[oiB s hape,<f>{<t>) + 2 B{(j))cos0 

- {sinOy 1 {P e - aP^)} (sinO^O), (7) 

where 

B shape, <f>{4>) = (mo/2 )M 2 Cl(N d — XX N d —yyj 1 

(8a) 

Pq = M [h x cosO coscj) + h y cosOsinf — h z sinO \, 

(8b) 

P(f> = M [h y coscj) — h x sincj)\ , (8c) 

hi= ^ ppSt G( °’ i} ( 8d > 

1/ At is proportional to the attempt frequency of the thermal 
field, At is the simulation time-step, and G( 0 ,i) is a Gaussian 
distribution with zero mean and unit variance (23) . 

III. Simulation Results 

We consider the magnetos trie tive layer to be made of poly¬ 
crystalline Terfenol-D, which has the following material prop¬ 
erties - magnetostrictive coefficient ((3/2)A s ): +90xl0 -5 , 
Young’s modulus (Y): 80 GPa, saturation magnetization 


(M s ): 8xl0 5 A/m, and Gilbert’s damping constant (a): 0.1 
(Refs. El, E3, d, G2). We use Terfenol-D (TbDyFe), 
a specially designed composite material rather than common 
magnetic materials (e.g., iron, nickel, or cobalt) since Terfenol- 
D has 30 times higher magnetostrictive coefficient in magni¬ 
tude 0 . Another good choice may be FeGa alloy, which has 
magnetostrictive coefficient +15x 10 -5 (28l . Since Terfenol-D 
has positive magnetostrictive coefficient, a compressive stress 
is required for magnetization to switch d . In our convention, 
a tensile stress is positive and a compressive stress is negative. 

The piezoelectric layer is considered to be made of lead- 
zirconate-titanate (PZT), which has a dielectric constant of 
1000 d. The PZT layer is assumed to be four times thicker 
than the magnetostrictive layer so that any strain generated 
in it is transferred almost completely to the magnetostrictive 
layer ©. The maximum strain on the PZT layer is considered 
to be 500 ppm l29l . t30l and it would require an electric field 
of 2.78 MV/m because d 3 i=L8xlO- 10 m/V for PZT ED. 
(We can also use relaxor ferroelectrics for the piezoelectric 
layer e.g., PMN-PT, PZN-PT that have high piezoelectric co¬ 
efficients and generate anisotropic strain, hence it can decrease 
the electric field further (32), (9), (T2|.) The corresponding 
stress in Terfenol-D is the product of the generated strain 
(500 x 10 -6 ) and the Young’s modulus (80 GPa). Hence, 
40 MPa is the maximum stress that can be generated in the 
Terfenol-D nanomagnet. For example, to generate 30 MPa 
stress, it would require voltages of 83.3 mV and 50 mV for 
nanomagnets with thicknesses 10 nm and 6 nm, respectively. 

We always ensure that the magnetostrictive nanomagnet has 
a single ferromagnetic domain EB, (33l for varied thickness 
and aspect ratio. The performance metrics switching delay 
and energy dissipation are determined by solving stochastic 
Landau-Lifshitz-Gilbert equation in the presence of room- 
temperature (300 K) thermal fluctuations following the pre¬ 
scription in Ref. @, f34l . We assume a 50 ps ramp for stress 
on the magnetostrictive nanomagnet 0. 

A. Lateral Dimensions and Demagnetization Factors 

Fig. |2(a)| plots the major axis a of the nanomagnet versus its 
thickness l with aspect ratio a/b as a parameter for an energy 
barrier height of 40 kT at room-temperature. The decrease in 
lateral dimensions with the increase in thickness are adjusted 
to satisfy the single-domain limit assumption ED . With these 
dimensions, the demagnetization factors Nd- yy and N d - Z z 
both increase with thickness and aspect ratio, however, N d - X x 
has the opposite trend ED- Fig- |2(b)| plots the difference 
Nd- yy — Nd-zz , which increases with both thickness and 
aspect ratio. Fig. |2(c)| plots N d - X x ~ Nd - yy , which decreases 
with the increase of both thickness and aspect ratio. 

B. Critical Stress and Success Rate of Switching 

Figs. |3(a)| and |3(b)| show the trend in lateral area and volume 
of the nanomagnet while varying its thickness with aspect 
ratio as a parameter. These plots can be derived from the 
Fig. |2(a)[ The decrease in lateral area bodes well with the 
requirement of scaling down the device size to enhance device 
density on a chip. However, the decrease in volume has to be 
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Barrier=40 kT Barrier=40 kT Barrier=40 kT 





Fig. 2. Major axis a of the nanomagnet, Nd_ yy — N^_ zz , and Ny_ xx — Nd- yy versus thickness l with aspect ratio a/b as a parameter to satisfy the 
single-domain limit for an energy barrier height of 40 kT at room-temperature, (a) Major axis a versus thickness l with aspect ratio a/b as a parameter. The 
major axis can be scaled down with the increase of both thickness l and aspect ratio a/b. (b) Nd_ yy — Ny_ zz versus thickness l with aspect ratio a/b 
as a parameter. Nd- yy — N^_ zz increases with both thickness l and aspect ratio a/b. (c) Nj-xx — Nd- yy versus thickness l with aspect ratio a/b as a 
parameter. Nd- xx — N c j_ yy decreases with the increase of both thickness l and aspect ratio a/b. 





Fig. 3. Elliptical cross-sectional area A of the nanomagnet, volume Q of the nanomagnet, and the minimum stress a min (to overcome the energy barrier 
height of 40 kT at room-temperature) versus thickness l with aspect ratio a/b as a parameter, (a) Elliptical cross sectional area A of the nanomagnet versus 
thickness l with aspect ratio a/b as a parameter. Area can be scaled down with the increase of both thickness l and aspect ratio a/b. (b) Volume Q of the 
nanomagnet versus thickness l with aspect ratio a/b as a parameter. Volume eventually scales down (note that area A scales down but thickness is scaled up) 
with the increase of both thickness l and aspect ratio a/b. (c) Minimum stress [crmin, see 0] versus th ickne ss l with aspect ratio a/b as a parameter. Note 
that crmin is proportional to N f j_ yy — N^_ zz and thus the curves follow the same trend as in the Fig. |2(b)| 


compensated by increasing the in-plane shape-anisotropy (oc 
Nd-yy — Nd-zz) to keep the energy barrier height same as 40 
kT (see ©). Thus, the curves in Fig. |2(b)| show the opposite 
trend to that of Fig. |3(b)| with the thickness and aspect ratio. 

Fig. |3(c)| shows the minimum stress cr m in required for 
stress-anisotropy to overcome the shape-anisotropy energy 
barrier of the nanomagnet. Since both the anisotropies are 
proportional to volume, crmin is independent of volume but 
it is proportional to Nd- yy — Nd- zz for a fixed barrier height 
of 40 kT (see 0 and Fig. |2(b)| ). 

Fig. [4] plots the successful switching rate versus stress for 
different thicknesses of the nanomagnet. A higher stress keeps 
the magnetization more out-of-plane (x-direction) and thus it 
is conducive to a higher success rate of switching f22l . For 
a fixed aspect ratio, a higher thickness of the nanomagnet 
corresponds to lower volume and thus lesser stress-anisotropy 
is produced at lower stress levels. This is the reason behind 
having lower success rate of switching with increasing thick¬ 
ness at the low stress levels. 


Barrier=40 kT a/b=l.l 



Fig. 4. Percentage of successful switching events (SR) as a function of stress 
(10-30 MPa) versus thickness l with aspect ratio a/b as a parameter for an 
energy barrier height of 40 kT. A moderately large number of simulations 
(10000) are performed in the presence of room-temperature (300 K) thermal 
fluctuations to generate each data point in these plots. At lower stress-levels 
(10-15 MPa), the number of successful switching events decreases with 
increasing thickness. 
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Barrier=40 kT Stress=30MPa 



Barrier=40 kT Stress=30MPa 



Fig. 5. Mean and standard deviation of switching delay (r mea n and r st d, 
respectively) versus thickness l with aspect ratio a/b as a parameter for an 
energy barrier height of 40 kT at room-temperature and stress 30 MPa. Each 
data point in these plots is calculated from 10000 simulations in the presence 
of room-temperature thermal fluctuations, (a) Thermal mean switching delay 
versus thickness l with aspect ratio a/b as a parameter. The mean switching 
delay increases with both thickness (for a fixed aspect ratio) and aspect 
ratio (for a fixed thickness), (b) Standard deviation in switching delay versus 
thickness l with aspect ratio a/b as a parameter. The standard deviation 
increases with both thickness (for a fixed aspect ratio) and aspect ratio (for a 
fixed thickness). 

C. Switching Delay and Energy Dissipation 

Figs. [5] and [6] plot the metrics switching delay and energy 
dissipation versus thickness with aspect ratio as a parameter 
for stress 30 MPa and an energy barrier height of 40 kT 
at room-temperature. In Fig. [3 we see that both the mean 
and standard deviation of switching delay become higher 
with increasing thickness. The reason behind the increase in 
switching delay is that it requires a minimum stress drain to 
overcome the energy barrier; drain increases with thickness 
(see Fig. |3(c)| ) and hence for a fixed stress of 30 MPa, 
the torque acting on the magnetization becomes lesser with 
higher thicknesses. The standard deviation in switching delay 
increases due to the decrease of nanomagnet’s volume with the 
increasing thickness (see Fig. |3(b)| ), which makes magnetiza¬ 
tion more prone to thermal fluctuations (see (l8dl0. The energy 
dissipation curves show an opposite trend because of the delay- 
energy trade-off. The “CV 2 ” energy dissipation as shown in 
Fig. |6(b)| decreases with increasing thickness as explained in 
the Section U 

IV. Discussions and Conclusions 

We have analyzed the area-delay-energy trade-offs of strain- 
coupled multiferroic composite devices. It turns out that it 


Barrier=40 kT Stress=30MPa 



Barrier=40 kT Stress=30MPa 



(b) 


Fig. 6. Total and “CV 2 ” energy dissipation ( E tota i and E, cv 2 ,, respec¬ 
tively) versus thickness l with aspect ratio a/b as a parameter for an energy 
barrier height of 40 kT at room-temperature and stress 30 MPa. Each data 
point in these plots is calculated from 10000 simulations in the presence 
of room-temperature thermal fluctuations, (a) Total energy dissipation versus 
thickness l with aspect ratio a/b as a parameter. Total energy dissipation 
increases with both thickness (for a fixed aspect ratio) and aspect ratio (for 
a fixed thickness), (b) “CV 2 ” energy dissipation versus thickness l with 
aspect ratio a/b as a parameter. “CV 2 ” energy dissipation decreases with 
the increase of both thickness (for a fixed aspect ratio) and aspect ratio (for 
a fixed thickness). 

is possible to increase the thickness of the magnetostrictive 
nanomagnet while reducing the lateral area of the nanomagnet 
keeping the single-domain limit assumption intact. However, 
the resultant volume of the nanomagnet eventually decreases 
with the decrease of lateral area. This weakens the induced 
stress-anisotropy in the nanomagnet for a material with fixed 
magnetostrictive coefficient and a certain stress value. The 
switching delay increases with the increasing thickness, while 
energy dissipation follows the opposite trend. An increase of 
stress can make the switching delay lower, however, there is 
a limit on stress that we can apply since it is dictated by 
the dielectric breakdown of the piezoelectric layer and the 
associated reliability issues. Thus, increasing magnetostrictive 
coefficient by superior material synthesis may be a solution 
to decrease switching delay as we reduce the lateral area 
(and increase the thickness) of a device. For example, us¬ 
ing a twice magnetostrictive coefficient of Terfenol-D (see 
Ref. 03, however, this is for single-crystal case rather than 
poly crystalline Terfenol-D) for thickness 10 nm, we can reduce 
the switching delay to the value same as for thickness 6 
nm, while achieving around two times reduction in energy 
dissipation. Both reducing lateral dimensions and increasing 
magnetostrictive coefficient are challenging in fabrication and 
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material synthesis procedures, which can increase the device 
density on a chip. 

Note that it may be possible to devise different switching 
methodologies with strain-mediated multiferroic composites 
having better performance metrics particularly reducing the 
lateral area of the devices further. However, since the strain 
generated is proportional to volume and with scaling volume 
reduces, it may be very useful to go beyond volumetric effect 
and employ surface-related phenomena for magnetoelectric 
coupling [36|, [37). 

With the growing experimental efforts, such devices may be 
a staple of modern non-volatile logic and memory systems for 
our future information processing paradigm and also could be 
deemed suitable for energy harvesting applications. 

Appendix A 

Fluctuations of Magnetization around Easy Axis 
due to Thermal Agitations 

Magnetization in a nanomagnet cannot move exactly from 
the easy axes (0 = 0°,180°) because the torque acting on 
magnetization become zero at these points [[22]. Fortunately, 
thermal fluctuations can deflect magnetization from these 
positions. Fig. [7 ] shows the distribution of polar angle 0 when 
magnetization is fluctuating around the easy axis 0 = 180° for 
an energy barrier height of 40 kT and for a nanomagnet with 
10 nm thickness and aspect ratio of 1.1. Similar Boltzmann 
distributions are achieved for different dimensions of the 
nanomagnet but for same barrier height 40 kT. The mean 
of this distribution is used as the magnetization’s initial polar 
angle to start with. 

Fig. [5] shows the distribution of azimuthal angle 4> when 
magnetization is fluctuating around the easy axis 6 = 180° for 
an energy barrier height of 40 kT. Two cases corresponding to 
two different thicknesses of the nanomagnet (6 nm and 10 nm) 
are considered. The distributions are Gaussian peaked with 
the peaks on the plane of the nanomagnet (0 = ±90°). We 
have assumed the initial value of cj) as 90°, however, assuming 
the other in-plane angle of magnetization produces similar 
magnetization dynamics due to symmetry of the underlying 
equations. Note that the distribution corresponding to the 10 
nm thickness one is less peaked since it is easier for the 
magnetization to deflect out-of-plane with higher thickness 
(Nd-xx is decreased at higher thickness, see Fig. |2(c)| ). 

Appendix B 

Results in Tabular Format 

Table H tabulates the simulation results presented in this 
paper. Different performance metrics e.g., successful switching 
rate, switching delay, and energy dissipation are given as a 
function of thickness with aspect ratio of the elliptical cross- 
section of the nanomagnet as a parameter. As we scale down 
the lateral dimensions of the nanomagnet, we can increase the 
thickness to some extent to keep the single-domain assumption 
valid. However, the volume still decreases with scaling and 
the volumetric strain decreases proportionately. Therefore the 
minimum stress cr m ^ n required to topple the energy barrier 
of the nanomagnet increases. Switching delay (both mean 


Distribution of e 



Fig. 7. Distribution of polar angle 9 when magnetization is fluctuating around 
the easy axis 9 = 180° for an energy barrier height of 40 kT. This plot 
is particularly for a nanomagnet with 10 nm thickness and aspect ratio of 
1.1. The most likely value of this distribution is 180° but the mean of this 
distribution is 175.55°. 


Distribution of <i> 



Distribution of <i> 



Fig. 8. Distributions of azimuthal angle 4> when magnetization is fluctuating 
around the easy axis 9 = 180° for an energy barrier height of 40 kT. (a) 
Distribution of azimuthal angle </> for a nanomagnet with 6 nm thickness and 
aspect ratio 1.1. (b) Distribution of azimuthal angle </> for a nanomagnet with 
10 nm thickness and aspect ratio 1.1. 

and standard deviation) increases and energy dissipation (both 
E tota i and 4 CV 2 ’ one) decreases as already discussed in the 
paper. 

For a given thickness, with the increase of aspect ratio, the 
anisotropy increases and it allows us to decrease the lateral 
area so that the energy barrier height can be kept constant. 
The volume decreases in turn and the minimum stress a m i n 
increases. Switching delay (both mean and standard deviation) 
increases and energy dissipation (both E tota i and ‘ CV 2 ’ one) 
decreases as already discussed in the paper. For thicknesses of 
9 and 10 nm, and aspect ratio of 1.2, we notice some switching 
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TABLE I 

Performance metrics e.g., success rate of switching (SR), mean switching delay (r mean ), standard deviation in switching delay 
(r std ), total energy dissipation CE toto i), and ‘CV 2, energy dissipation^^, for a Terfenol-D/PZT STRAIN-MEDIATED MULTIFERROIC 

COMPOSITE DEVICE WITH RESPECT TO THICKNESS l WITH THE ASPECT RATIO ( a/b ) OF THE ELLIPTICAL CROSS-SECTION OF THE NANOMAGNET AS A 
PARAMETER. THE STRESS IS 30 MPA AND THE ENERGY BARRIER HEIGHT IS 40 kT . THE DERIVED QUANTITIES THAT ARE GIVEN ARE 
N d _y Z = N d _yy — N d _ zz , N d _ X y = N d _ xx — N d _yy, Area (A), VOLUME (Q), AND MINIMUM STRESS TO OVERCOME THE ENERGY BARRIER OF 
THE NANOMAGNET (7min. THE PERFORMANCE METRICS ARE CALCULATED FROM 10000 SIMULATIONS. 


1 

(nm) 

a/b 

a 

(nm) 

N d -y Z 

N d — xy 

A 

(nm 1 2 ) 

Q 

(nm 3 4 5 6 7 8 9 10 ) 

®min 

(MPa) 

SR 

(%) 

Tmean 

(ns) 

T std 

(ns) 

Rtotal 

(kT) 

^‘CE 2 / 

(kT) 

6 

1.10 

100.0 

0.0095 

0.7754 

7139.98 

42839.90 

4.25 

100.00 

0.35 

0.042 

328.45 

49.85 


1.15 

76.5 

0.0170 

0.7206 

3996.82 

23980.94 

7.60 

100.00 

0.38 

0.063 

169.95 

15.66 


1.20 

64.0 

0.0253 

0.6756 

2680.83 

16084.95 

11.32 

100.00 

0.42 

0.080 

113.40 

7.05 

7 

1.10 

82.0 

0.0121 

0.7112 

4800.92 

33606.47 

5.41 

100.00 

0.36 

0.051 

237.05 

22.59 


1.15 

63.0 

0.0214 

0.6449 

2710.65 

18974.54 

9.55 

100.00 

0.42 

0.076 

130.51 

7.21 


1.20 

53.0 

0.0315 

0.5921 

1838.49 

12869.40 

14.05 

100.00 

0.49 

0.094 

92.11 

3.32 

8 

1.10 

69.5 

0.0147 

0.6456 

3448.79 

27590.32 

6.57 

100.00 

0.39 

0.065 

186.85 

11.67 


1.15 

54.0 

0.0255 

0.5712 

1991.50 

15931.97 

11.41 

100.00 

0.47 

0.087 

110.15 

3.89 


1.20 

45.5 

0.0372 

0.5112 

1354.98 

10839.80 

16.63 

100.00 

0.57 

0.114 

81.65 

1.80 

9 

1.10 

60.5 

0.0172 

0.5806 

2613.41 

23520.71 

7.70 

100.00 

0.43 

0.076 

157.00 

6.71 


1.15 

47.0 

0.0296 

0.4968 

1508.65 

13577.83 

13.23 

100.00 

0.54 

0.098 

96.55 

2.24 


1.20 

40.0 

0.0426 

0.4335 

1047.20 

9424.78 

19.04 

99.99 

0.68 

0.139 

76.61 

1.08 

10 

1.10 

53.5 

0.0197 

0.5159 

2043.64 

20436.42 

8.79 

100.00 

0.48 

0.081 

136.85 

4.10 


1.15 

42.0 

0.0333 

0.4275 

1204.73 

12047.32 

14.88 

100.00 

0.62 

0.117 

89.59 

1.43 


1.20 

36.0 

0.0474 

0.3620 

848.23 

8482.30 

21.18 

99.91 

0.81 

0.168 

75.44 

0.71 


failures among the 10000 simulations due the increase in the 
minimum stress a min- Note that although the stress on the 
nanomagnet 30 MPa is higher than crmin, sometimes thermal 
fluctuations can scuttle the magnetization causing switching 
failures. 
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